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Abstract 

A simple statistical mechanical approach for studying multilayer adsorption of in- 
teracting polyatomic adsorbates (/c-mers) has been presented. The new theoretical 
framework has been developed on a generalization in the spirit of the lattice-gas 
model and the classical Bragg- Williams (BWA) and quasi-chemical (QCA) approx- 
imations. The derivation of the equilibrium equations allows the extension of the 
well-known Brunauer-Emmet- Teller (BET) isotherm to more complex systems. The 
formalism reproduces the classical theory for monomers, leads to the exact statisti- 
cal thermodynamics of interacting k-mers adsorbed in one dimension, and provides 
a close approximation for two-dimensional systems accounting multisite occupancy 
and lateral interactions in the first layer. Comparisons between analytical data and 
Monte Carlo simulations were performed in order to test the validity of the theo- 
retical model. The study showed that: (i) the resulting thermodynamic description 
obtained from QCA is significantly better than that obtained from BWA and still 
mathematically handable; (ii) for non-interacting fc-mers, the BET equation leads 
to an underestimate of the true monolayer volume; (Hi) attractive lateral interac- 
tions compensate the effect of the multisite occupancy and the monolayer volume 
predicted by BET equation agrees very well with the corresponding true value; 
and (iv) repulsive couplings between the admolecules hamper the formation of the 
monolayer and the BET results are not good (even worse than those obtained in 
the non-interacting case). 

Key words: Equilibrium thermodynamics and statistical mechanics, Surface 
thermodynamics, Adsorption isotherms, Monte Carlo simulations 
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1 Introduction 



Adsorption on solid surfaces is a complex phenomenon [1,2,3,4,5], which im- 
plies a series of questions about the nature of the forces binding foreign 
molecules to a surface and about the thermodynamic behavior of the sys- 
tem. Although the problem of the forces is still far from being elucidated, a 
large amount of interest has been devoted to the study of statistical proper- 
ties, within the frame of simplified models for both monolayer and multilayer 
adsorption. In 1918 Langmuir derived the monolayer adsorption isotherm ki- 
netically for gas molecules adsorbed on the homogeneous surface of adsorbents 
without attractions among the adsorbed molecules [6]. After that Brunauer, 
Emmett and Teller [7], based on a model of localized adsorption, developed 
the most important theory of multilayer adsorption. Their equation, the BET 
isotherm, was the first and the most useful, covering the complete range of 
pressures up to p , the saturation pressure. Later, Hill [1] derived the BET 
isotherm statistically on a group of homogeneous adsorption sites for the mul- 
tilayer adsorption since it was derived kinetically by Brunauer, Emmett and 
Teller. It is found to be in good agreement with some experimental data for 
relative pressures less than about 0.5 [8]. 

But the theoretical BET isotherm deals with the common assumption that 
each ad-molecule occupies one adsorption site of the surface. However, most 
adsorbates involved in experiments are polyatomic; hence, the theoretical de- 
scription of their thermodynamic properties is a topic of much interest in 
adsorption theoryt!]- Leading contributions to this subject, generically called 
multisite-occupancy adsorption, have been presented by Flory [9], Huggins 
[10], Guggenheim [11], DiMarzio [12], Nitta et al. [13] and Rudzinski et al. [14] 
through approximate treatments of monolayer adsorption on homogeneous 
and heterogeneous surfaces. More recently, Aranovich and Donohue [15,16] 
derived a multilayer adsorption isotherm, which is not limited by the func- 
tional form of the monolayer adsorption isotherm and should be capable to 
include multisite occupancy (with an adequate choice of a fitting parameter). 
On the other hand, the closed exact solution for the multilayer adsorption 
isotherm of dimers, along with the basis for calculating adsorption thermody- 
namics of homonuclear polyatomic molecules (fc-mers) on one-(lD) and two- 
dimensional (2D) substrates, have been recently presented [17,18]. This rig- 
orous thermodynamic study demonstrated that the entropic contribution of 
non-spherical adsorbates is significant in the multilayer regime when compared 
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2 Even for simple gases such as oxygen, nitrogen and carbon monoxide, which 
basically are not altered in their molecular dimensions under physical adsorption, 
the adsorption energy depends in general on the orientation of the molecule in the 
adsorbed state. 
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with monoatomic adsorption. Thus, the determinations of surface areas and 
adsorption energies from polyatomic adsorbate adsorption may be severely 
misestimated, if this polyatomic character is not properly incorporated in the 
thermodynamic functions from which experiments are interpreted. 

There is another important physical fact which has not been sufficiently stud- 
ied; namely, the effect of the lateral interactions between the ad-molecules 
in presence of multilayer adsorption and multisite occupancy. Thus, most of 
the theoretical work dealing with adsorption of interacting polyatomics has 
been based on models of monolayer adsorption. In this context, the aim of 
the present work is to extend the treatment of Refs. [17,18] to include lat- 
eral interactions in the adsorbate. Here, we introduce nearest-neighbor in- 
teractions between the molecules adsorbed in the first layeiGO, by following 
the configuration-counting procedure of the Bragg- Williams approach and the 
quasi-chemical approximation. In addition, Monte Carlo (MC) simulations are 
performed in order to test the validity of the theoretical model. The new the- 
oretical scheme allows us: (1) to reproduce the classical theory for monomers 
[1,5]; (2) to develop a closed exact expression for the multilayer adsorption 
isotherm of interacting /c-mers on ID chains; (3) to obtain an accurate ap- 
proximation for multilayer adsorption on 2D substrates accounting multisite 
occupancy and lateral interactions; and (4) to provide a simple model from 
which experiments may be reinterpreted. 

The present work is organized as follows. In Section 2, the theoretical formal- 
ism along with the basis of the MC method are presented. In Section 3, the 
results of the theoretical model are shown and discussed by comparing with 
MC simulations. Finally, conclusions are drawn in Section 4. 



2 Basic Definitions: Adsorption Model and Monte Carlo Simula- 
tion 

2.1 Model 

In this section we present the lattice-gas model for the adsorption of particles 
with multisite occupancy in the multilayer regime. The adsorbent is a homo- 
geneous lattice of sites with coordination number (. The adsorbate is assumed 

3 As it is well-known, BET equation can be applied at coverage not greatly ex- 
ceeding (statistically) monolayer coverage. Thus, although the contribution from 
the secondary adsorption can already be significant, the density of the molecules in 
the second and higher adlayers is expected to be much lower than that in the first 
adsorbed layer. Therefore, it seems to be satisfactorily enough to take into account 
only the interactions between the primarily adsorbed molecules [1]. 
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as linear molecules having /c-identical units (/c-mers) each of which occupies an 
adsorption site. Furthermore, i) a k-mer can adsorb exactly onto an already 
adsorbed one; ii) attractive and repulsive lateral interactions are considered 
in the first layer and horizontal interactions are ignored in higher layers; iii) 
the adsorption heat in all layers, except the first one, equals the molar heat 
of condensation of the adsorbate in bulk liquid phase. Thus, c = qi/qi = qi/q 
with qi = q (i = 2,. ..,00) denotes the ratio between the single-molecule 
partition functions in the first and higher layers [17]. The fact that fc-mers can 
arrange in the first layer leaving sequences of / empty sites with I < k, where 
no further adsorption of a k— mer can occur in such a configuration, makes the 
calculation of entropy much elaborated than the one for monomer adsorption. 

To describe a system of N A;-mers adsorbed on M sites at a given temperature 
T, let us introduce the occupation variable which can take the values s$ = 
or 1, if the site % is empty or occupied by a /c-mer unit, respectively. Under 
these conditions, the Hamiltonian of the system is given by 

M 

H = dSi + ke(N - N ± ) + ^ w$iSj - N ± (k - l)w, (1) 

where e\ (e) represents the adsorption energy of a A;-mer unit on the first layer 
(higher layers); Aq = s i/k is the number of /c-mers adsorbed on the first 
layer; w is the lateral interaction energy between two nearest-neighbor (NN) 
units belonging to different k-mers adsorbed in the first layer (we use w > 
for repulsive and w < for attractive interactions) and represents pairs 
of NN sites. The term Ni(k — l)w is subtracted in equation (1) since the 
summation over all the pairs of NN sites overestimates the total energy by 
including N ± (k — 1) bonds belonging to the Ni adsorbed k-mers. 



2.2 Theory 

From a theoretical point of view, when intermolecular forces are introduced 
(in our case, NN interactions in the first layer), an extra term in the partition 
function for interaction energy is required. With this extra term, only partition 
functions for the whole system can be written. Ising [19] gave an exact solution 
to the ID monolayer in 1925. All other cases are expressed in terms of series 
solution [1,20], except for the special case of 2D monolayers at half-coverage, 
which was exactly solved by Onsager [21] in 1944. Close approximate solutions 
in dimensions higher than one can be obtained, and the two most important of 
these are the Bragg- Williams approximation (BWA) [1] and the quasi-chemical 
approximation (QCA) [1,22]. These leading models have played a central role 
in the study of adsorption systems in presence of lateral interactions between 
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the adatoms. Next, we apply BWA and QCA to study multilayer adsorption 
of interacting fc-mers. 



The BWA is the simplest mean-field treatment for interacting adsorbed parti- 
cles, even in the case of multilayer adsorption and multisite occupancy. For a 
lattice having M adsorption sites, the maximum number of columns that can 
be grown up onto it is n max = M/k. If an infinite number of layers is allowed 
to develop on the surface, the grand partition function, S, of the adlayer in 
equilibrium with a gas phase at chemical potential \i and temperature T, is 
given by 



~ = CMn, M, C) exp \-(3E k (n, M) 



(2) 



n=0 



where £ is the grand partition function of a single column of /c-mers having 
at least one k-mei in the first layer; Q k (n, M,() is the total number of dis- 
tinguishable configurations of n columns on M sites with connectivity ( and 
E k (n, M) is the mean total energy of the system assuming that the n columns 
are randomly distributed over the lattice. 

Then, 

£ = T i q 1 q i - 1 \ i = cfVA* = -^L = _£^_ (3) 



where A = exp(/x//c B T) is the fugacity and fc B is the Boltzmann constant. 
In addition, it is possible to demonstrate that x = \q = p/p Q is the relative 
pressure [5,17]. 

M, () can be approximated considering that the columns are distributed 
completely at random on the lattice and assuming the arguments given by 
different authors [23,24,25] to relate the configurational factor Qk{n, M,Q for 
any (, with the same quantity in the ID case (£ = 2). Thus 

n k (n,M,() = v(Ck) n n k (n,M,2), (4) 



where 77(^5 k) is, in general, a function of the connectivity and the size of 
the molecules. In the particular case of rigid straight /c-mers it follows that 
T]((, k) = C/2. In addition, f2fc(n, M, 2) can be readily calculated [25] giving 



On the other hand, E k (n, M) needs to be calculated. For this purpose, let us 
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consider the number of NN of a k-mer (column) adsorbed in the first layer, 
* = [2(C-l) + (*-2)(C-2)], (6) 



where the first term in the RHS of Eq. (6) is the number of NN connected to 
both extrems of the k-mer and the second term is due to the NN connected 
to the k — 2 inner units of the molecule. The probability that one of the NN is 
filled by molecules is equal to kn/M (random distribution of molecules among 
sites). Then, the mean number of pairs of filled sites is given by 



E k (n,M) = -n [2(( - 1) + (k 

zkn 2 
= 2M ' 



2)(C-2)] 



' kn 
~M 



(7) 



Then, the grand partition function can be written as, 

„ [M- (k- l)n}\ ( (3wzkn 2 \ 



When w — 0, the summation in Eq. (8) can be performed very easily. In the 
more general case w^O, the summation cannot be done in an easy way, but 
we can apply the method of the maximum term. Then, the sum in Eq. (8) 
is replaced by its maximum term, found from the condition = 0, being 

t(n, M, T) = [ V ^ ^fe^ exp (-^) . Thus, 

ln?7 f - P W ^ Ul — (k — l)ln[M- (k - l)m] -lnni + k\n(M - kn x ) = 0, 

where n\ is the value of n giving the maximum term in the sum in Eq. (8). 

Introducing n\ in Eq. (8) and applying Stirling's approximation, In S is given 
by 



In S = m In 77 e - + [M-(k- l)ni] ln[M — (k — l)m] 

-nilnm - (M- fcni)ln(M- fcni). (10) 



The average number of the molecules in the adsorption system N is 

N = k B T[ —- 

V J TM 
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Then, 



s 9fi ) TM 



-ri x lnni + kn[ ln(M — krii), 



where n[ = ^ and £' has the following explicit form: 



, = d^ = 1 qiexp(fi/k B T) 



dpi k B T[l- q ex.p(n/k B T)} 



2 ' 



(12) 



(13) 



Inserting the condition Eq. (9) into Eq. (12), we obtain 



'dlnE\ _ f 



(14) 



Considering now Eqs. (9), (11), (13) and (14), the adsorption isotherm equa- 
tion can be obtained. In the case of adsorbed monomers (k — l,rj — 1), Eq. 
(9) can be written as 



ln£ - - lnm + ln(M - m) = 0, 



(15) 



and 



0i 



rii £exp(— (5wz6\) 
M ~ l + ^exp(-(3wz6 i y 



where 6\ = kn\jM is the monolayer coverage, being 9 
coverage. Then, from Eqs. (3), (11), (13) and (16) 



(16) 

kN/M the total 



T,M 



1 - qexp(jj,/k B T)_ 

qi exp(fj,/k B T — wz9i/k B T) 



1 — qexp(—jj/k B T) + qi exp(jj/k B T — wz9\jk B T) 



(17) 



which can be easily recognized as the classical BET equation if we write it in 
the form 
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8 = ^—,i — 7, (18) 

(1 — x) (1 — X + c*x) 

where c* = (q 1 /q)exp(—wz9i/k B T). 

For the case of dimers (k — 2), Eq. (9) reduces to 

In 77 £ - - lnm - ln(M - m) + ln(M - 2m) 2 = 0, (19) 



and 



e 1 = ^± = i 1 (20) 

M [l + 4r/£exp(- / 3u^ 1 )] 1/2 



Using Eqs. (3), (11), (13) and (20) we obtain 



s = m=2k B T fa ^ 



M \ dpi 



T,M 



f l~ 1 - rTTTTs l (21) 



1 -qexp(/i/k B T) \ [i + 4r/ eexp(-/?w^i)] 1/2 
and, in terms of c* = r/(qi/q) exp(— wzOi/ksT) and x, 



= . 



1 — rr 



Hi/2' 



'1 + 4c*x — x) 



(22) 



Eq. (22) is similar to the recently reported multilayer isotherm for non-interacting 
dimers [18]. In this case, taking into account the lateral interactions between 
the primarily adsorbed molecules adds to the adsorption energy the term 
(3wz0i, which represents the potential of the average force acting on an ad- 
molecule in the first adsorbed layer from its NN in the first layer. 

For k > 2 the explicit expression of the adsorption isotherm cannot be ob- 
tained in a easy way. However, the calculations for large molecules can be 
easily done through a standard computing procedure; in our case, we used 
Maple software. 

An alternative method to calculate the multilayer adsorption isotherm was 
recently reported in Ref.[18]. The theoretical procedure can be described as 
follows: 

) By using 6\ as a parameter (0 < 6i < 1), the relative pressure is obtained 
by using the condition 



1 
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X 



1 

1 + cAi 



i > 



(23) 



where Ai is the monolayer fugacity. This calculation requires the knowledge 
of an analytical expression for the monolayer adsorption isotherm, Ai(#i). 
2) The values of 9\ and x are introduced in 



e 



9i 



1-x' 

and the total coverage is obtained. 



(24) 



The equivalence between both methodologies can be easily understood. In 
fact, Eq. (23) can be obtained from Eq. (3) and the maximum term condition 
Eq. (9). On the other hand, from Eqs. (11), (13) and (14) and after some 
algebra, the total coverage can be written in terms of the monolayer coverage, 
and Eq. (24) is recovered. As an example, in the following we show the use of 
the method in Ref.[18] to calculate the adsorption isotherms in Eqs. (18) and 
(22). 



We start from the equation 



Ai = 



0i 



fc-i 



exp (/3zw9i), 



(25) 



which represents the BWA isotherm of interacting /c-mers adsorbed at mono- 
layer [26,27]. 

Substituting Eq. (25) into Eq. (23), one obtains the following expression for 
the relative pressure, 



V 



9i 



1 



(fe-i). 



k-l 



9i exp ((3zw6i) 



Po C77(c, k)k (i - + e 1 



1 _ (fc- 1 ) / 



k-l 



9\ exp ((3zw9i] 



(26) 



Eqs. (24) and (26) represent the mean-field solution describing the adsorption 
of interacting /c-mers at multilayer regime on a homogeneous surface. In the 
case of monomer [dimer] adsorption, Eqs. (24) and (26) reduce to Eq. (18) 
[(22)]. 

We now turn to the QCA, which is significantly better than the BWA. The 
important assumption in this method is that pairs of NN sites are treated as 
if they were independent of each other (this assumption is, of course, not true, 
because the pairs overlap [1]). 
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In order to apply the scheme in Ref. [18], we start with the monolayer ad- 
sorption isotherm of interacting /c-mers adsorbed on a lattice of connectivity 
C obtained from the formalism of QCA [27], 



Ai 



0i exp (pwz/2) 



lfc-1 



2frj 

2fc 



— a 



z/2 



- (k - 1)9,) 



fc-1 



i(l-fli) -a 



^ (^) ; 



,(27) 



where a is 



a = 



2k 



(¥)»i + & 



2 I fc 



n 2 



01 



^40i(l-0i) 



1/2 



(28) 



(29) 



and 



A = 1 - exp(-^ty). 



(30) 



Replacing Eq. (27) into Eq. (23), we obtain 

,2(fc-l) r 



p_ 

3>o. 



1 + 



ckrj((, k) 



2 



-(k-m) 



k-l r 



Id 



a 



kC,/2 



0i exp ((3wz/2)(l - 0i) fc (C-!) [k — (k — l)0i 



fc-i 



201 

2fc 



2/2 



,(31) 



a 



Eqs. (24) and (31) represent the solution describing the multilayer adsorption 
of interacting fc-mers on homogeneous surfaces in the framework of the QCA. 



2.3 Monte Carlo Simulation of Adsorption in the Grand Canonical Ensemble 



The adsorption process is simulated through a grand canonical ensemble Monte 
Carlo (GCEMC) method[18]. 

For a given value of the temperature T and chemical potential //, an initial 
configuration with iV /c-mers adsorbed at random positions (on kN sites) 
is generated. Then, an adsorption-desorption process is started, where each 
elementary step is attempted with a probability given by the Metropolis [28] 
rule: 

W = mm {1, exp [-13 (AH - /2 AN)}} , (32) 
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where AH and AN represent the difference between the Hamiltonians and the 
variation in the number of particles, respectively, when the system changes 
from an initial state to a final state. In the process there are four elemen- 
tary ways to perform a change of the system state, namely, adsorbing one 
molecule onto the surface, desorbing one molecule from the surface, adsorbing 
one molecule in the bulk liquid phase and desorbing one molecule from the 
bulk liquid phase. In all cases, AN = ±1. 

The algorithm to carry out one MC step (MCS), is the following : 

1) Set the value of the chemical potential \i and the temperature T. 

2) Set an initial state by adsorbing iV molecules in the system. Each /c-mer 
can adsorb in two different ways: %) on a linear array of (k) empty sites on 
the surface or %%) exactly onto an already adsorbed /c-mer. 

3) Introduce an array, denoted as A, storing the coordinates of n e entities, 
being n e , 

n e = number of available adsorbed k— mers for desorption (n<}) 



where n a is the sum of two terms: i) the number of /c-uples of empty sites 



on the surface and ii) the number of columns of adsorbed fc-mera 4 1. 
4) Choose randomly one of the n e entities, and generate a random number 



4.1) if the selected entity is a /c-uple of empty sites on the surface then adsorb 
a /c-mer if £ < W™£ , being W™£ the transition probability of adsorbing 
one molecule onto the surface. 

4.2) if the selected entity is a /c-uple of empty sites on the top of a column of 
height i, then adsorb a new /c-mer in the i + 1 layer if £ < Wf™[ l s k , being 
W%$j, the transition probability of adsorbing one molecule in the bulk 
liquid phase. 

4.3) if the selected entity is a /c-mer on the surface then desorb the /c-mer 
if £ < Wfo s , , being W^l the transition probability of desorbing one 
molecule from the surface. 

4.4) if the selected entity is a /c-mer on the top of a column then desorb the 
/c-mer if £ < W^ k , being W^ k the transition probability of desorbing 
one molecule from the bulk liquid phase. 

5) If an adsorption (desorption) is accepted in 4), then, the array A is updated. 

6) Repeat from step 4) M times. 

In the present case, the equilibrium state could be well reproduced after dis- 
carding the first m ~ 10 6 MCS. Then, averages were taken over ml ~ 10 6 MCS 



4 Note that the top of each column is an available /c-uple for the adsorption of one 
/c-mer. 



+ number of available k— uples for adsorption (n a ), 



(33) 




£e [0,1] 



n 



successive configurations. The total coverage was obtained as simple averages, 



where (N) is the mean number of adsorbed particles, and (...) means the time 
average over the MC simulation runs. 

The computational simulations have been developed for ID chains of 10 4 sites, 
and square L x L lattices, with L = 100, and periodic boundary conditions. 
With this lattice sizes we verified that finite-size effects are negligible. 



3 Results 



In the present section, we will analyze the main characteristics of the ther- 
modynamic functions given in Subsection 2.2, in comparison with simulation 
results for a lattice-gas of interacting k-mers on ID and 2D substrates. 



3.1 Exact solution in the ID case 



In Figs. 1-3 we address the comparison between the analytical adsorption 
isotherms for ID substrates and MC simulations. Different values of the pa- 
rameter c, the lateral interactions and the A;-mer sizes have been considered. 

We start analyzing the case of k — 1, c = 1 and different values of w/UbT = 
-2, -1, 0, 1 and 2 [see Fig. 1 (a)]. The case w = (standard BET model) has 
been widely discussed in the literature (see, for instance, Ref. [8]) and it has 
been shown that a shape of a Type II isotherm is obtained so long as c exceeds 
2. In the case of this figure, c = 1 and the curve for non-interacting particles 
has the general shape of a Type III isotherm. For repulsive couplings, the 
interactions do not favor the adsorption on the first layer and the isotherms 
shift to higher values of pressure. On the other hand, attractive lateral inter- 
actions facilitate the formation of the monolayer. Consequently, the isotherms 
shift to lower values of p/p and their slope increases as the ratio \w\/kBT 
increases. With respect to the shape of the curves, there exists a range of 
w/k B T where the isotherms keep the shape of a Type III isotherm (in this 
case, —1 < w/k B T < 1). However, a knee appears in the isotherms (the curves 
adopt the shape of a Type II isotherm) as the ratio \w\/kBT increases. The 
shape of the knee depends on the value of \w\/kBT, becoming sharper as the 
value of w/ksT becomes more negative. This point will be illustrated more 
clearly in Fig. 5. 
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The effect of the k-mer size on the adsorption isotherm can be understood 
by analyzing Figs. 1 (b) and 1 (c), where the study in Fig. 1 (a) is repeated 
for k = 2 and k — 4, respectively. Two main conclusions can be drawn from 
the figures. Namely, 1) the difference between the curves corresponding to 
different values of w/ksT and 2) the range of w/k B T where the isotherms do 
not develop an inflection point diminish as k is increased. 

Figs. 2 and 3 show the effect of the parameter c on the adsorption isotherms. 
As can be observed, all curves exhibit a pronounced knee as the parameter c 
is increased. This effect can be better visualized in the insets of Figs. 2 and 
3, where a zoom of the region of low pressure is presented. In the case of 
attractive interactions, the knee appears around 9 = 1 and can be associated 
to the formation of the monolayer. In the case of repulsive interactions, /c-mers 
avoiding configurations with NN heads arrange in a structure of alternating 
particles separated by an empty site. Thus, for a given value of k and strong 
repulsive interactions, a marked knee is found at 9 = kjik + 1). 

To complete the discussion of Figs. 1-3, we evaluate the reaches and limitations 
of the two theoretical approximations studied. QCA leads to exact solution in 
ID systems. Consequently, MC simulations in the grand canonical ensemble 
(symbols) fully agree with the predictions from QCA (solid lines), which rein- 
forces the robustness of the two methodologies employed here. With respect 
to BWA, two different behaviors are observed : {%) for small values of |u>|//cbT, 
the theoretical curves show a good agreement with the simulation data and 
{%%) for \w\/k B T > 2, appreciable differences are observed between BWA and 
MC results. For strong attractive couplings [see, for instance, the case corre- 
sponding to w/ksT = —2 in Fig. 1 (a)], a characteristic van der Waals loop 
is observed in the adsorption isotherm and BWA incorrectly predicts a phase 
transition for ( = 2. The shape of the isotherms is fairly independent of the 
size of the molecules. However, the disagreement between the BWA curves 
and the exact results turns out to be significantly large for larger k-mer sizes 
[see, for instance, Fig. 1 (c)]. 

The analysis of the curves in Figs. 1-3 indicates that the appearance of a 
inflection point in the curves depends on k, c and w/ksT. This will be studied 
in detail in the following. The point of inflection can be obtained in three 
steps: (1) differentiating twice the adsorption isotherm equation to obtain 
d 2 9/dY 2 (being Y = p/p for the sake of simplicity); (2) equating the resulting 
expression to zero and solving for Y gives Y F , the value of p/p at the point 
of inflection; and (3) inserting Y F in the adsorption isotherm equation gives 
9 F) the value of 9 at the point of inflection. 

The location of the point of inflection (X F = 9 F , Y F ) is plotted in Fig. 4 for 
different values of k and w/ksT. The information is organized as follows: (i) 
as in Figs. 1-3, parts (a), (b) and (c) correspond to k = 1, k = 2 and k = 4, 
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respectively; (ii) the curves in (a), (b) and (c) were obtained for different 
values of w/ksT (as indicated in the caption of each figure); and (Hi) each 
point on a given curve corresponds to a determined value of c. 

In order to understand the basic phenomenology, we consider in the first place 
the case corresponding to w/k B T = (highlighted curve). Clearly, the value 
of 9 at the point of inflection may deviate considerably from unity. However, 
there exist a certain value of c = c m , where the point of inflection coincides 
with the point corresponding to the monolayer capacity. Fig. 5 shows the 
values of c m , obtained numerically, as a function of w/ksT for four k-mer 
sizes (k = 1,2,4,10). Two regimes can be clearly differentiated according 
to the sign of the lateral interactions. For attractive interactions, c m is not 
defined in all the range of w/ksT. Thus, for each k-mer size, there exists a 
limit value of the lateral interaction, w m j n /fc^T, below of which the coordinate 
Xp characterizing the inflection point in the adsorption isotherm is larger than 
one. In other words, the Xp — Yp diagrams corresponding to values of w/ksT 
below Wmin/kBT do not cross the line corresponding to Xp = 1 (dashed line 
in Fig. 4). The values of w min /k B T for k — 1, 2, 4, 10 are collected in Table I. 
Finally, c m increases monotonically as the interaction energy is increased in the 
range w min /k B T < w/ksT < 0. On the other hand, c m shows an exponential 
dependence [exp(nw /ksT)] in the range w/ksT > 0, where k is a parameter 
depending on the k-mer size. The different values of k, obtained from the slope 
of lnc m vs. w /ksT are reported in Table I. 

For values of c between c m and infinity the adsorption at the point of inflection 
exceeds the monolayer capacity; for values of c below c m the two quantities 
deviate more and more and for a limit value of c = c n , the point of inflection 
disappears^. In the low-coverage regime (9 — ► 0), c n vs. w/k B T can be cal- 
culated analytically (see Appendix for further discussion). The result of this 
calculation is presented in Fig. 6 for different values of k. Solid lines represent 
theoretical data from Eq. (73) and symbols correspond to values of c n obtained 
numerically. As can be visualized from the figure, the dashed line separates 
two well differentiate regions. At right of the dashed line, the coordinates of 
the inflection point corresponding to the limit value c n are Xp = and Yp = 0. 
Then, the assumption of 9 — > in Appendix is valid and the symbols coincide 
with the solid line. At left of the dashed line, the point of inflection disappears 
for Xp > and the solid line is not defined. The value of w/ksT correspond- 
ing to the dashed line was obtained from the condition 2k + 1 — 2e~ M,//fcsT = 
[see Eq. (73)]. From the point of view of the X F — Y F diagrams, the previous 
condition separates diagrams defined in the origin (Xp = 0,Yp = 0) from 
those where the inflection point disappears for (Xp > 0,Yp > 0) (see solid 
circles in Fig. 4). 



5 For c > c n , the isotherm is of Type II and when c is less than c n the isotherm is 
of Type III and discussion of the point of inflection is meaningless. 
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In the next we will refer to one of the main applications of BET model, which 
consists in taking an experimental isotherm in the low-pressure region and 
fitting values of the monolayer volume and the parameter c, from the linearized 
form of the BET equation, 

p/p. _ i + tzH p/p .. (35) 



[v (1 - p/p 



The plot of (p/p )/[v(l — p/Po)} v s p/p should therefore be a straight line 
with slope s = (c — l)/cv m and intercept i = l/cv m . Solution of these two 
simultaneous equations gives v m and c: 

1 s 

v m = and c = — hi. (36) 

s + 1 i 



In this context, it is of interest to study the behavior of multilayer isotherms 
of fc-mers (with k > 2) in the low-pressure regional in comparison with BET 
isotherm. For this purpose, we will analyze, by using the standard BET formal- 
ism, exact theoretical isotherms in the ID case. As an example, Fig. 7 shows 
the results obtained for k = 2, c = 10, w/k B T = —2, —1, 0, 1, 2 and pressures 
ranging from p/p Q = up to p/p = 0.30. Symbols represent theoretical data 
from QCA (exact results) and lines correspond to linearized forms of the BET 
equation. A linear function is only obtained if k = 1 and w/ksT = (see 
inset). The nonlinear behavior of interacting /c-mers isotherms at low pres- 
sures, which is a distinctive characteristic of many experimental isotherms, is 
showing that the polyatomic character of the adsorbate and the lateral inter- 
actions must be taken into account. The significant differences observed as k 
and w/ksT are varied indicate that the analysis of experimental isotherms 
of interacting larger molecules by means of the BET isotherm would lead to 
values of the parameters c and v m appreciably different from the real ones. 

In order to measure the differences mentioned above, the analysis of Fig. 7 was 
repeated in Fig. 8 for c = 10 [part (a)], c = 100 [part (b)] and different values 
of k and w/ksT. The results obtained for the monolayer volume are shown 
in the figure, where v m and vbet represent the real monolayer volume and 
the corresponding value obtained from the BET fitting, respectively. We start 
analyzing the case of w/ksT = (see inset). In this condition, the difference 
between v m and vbet (^bet/ ^m) increases (decreases) initially as the /c-mer 
size is increased and remains almost constant for larger values of k. As is shown 
in the figure, these differences diminish as the parameter c increases. Now, it 

6 Although in each particular case it is possible to find an optimum range of relative 
pressures, for practical purposes, we have chosen to set this range from 0.05 to 0.25. 
Nevertheless, by choosing other ranges (for example, between 0.05 and 0.35) we 
obtain similar results. 
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is interesting to analyze the effect of the lateral interactions. As was discussed 
above, attractive lateral interactions favor the formation of the monolayer 
and, consequently, compensate the effect of the multisite occupancy. Thus, 
for a given value of k, vb et / % tends to one as \w\/ksT is increased. On 
the other hand, repulsive interactions do not facilitate the formation of the 
monolayer, increasing the differences between v m and v B et- For w/ksT » 1, 
a marked knee is found at 9 — k/{k + 1) and vbet is close to this value of 
coverage. As the fc-mer size is increased, 9 — k/(k + l) — > 1 and, consequently, 
VBEr/vm 1- This phenomenon can be clearly visualized by observing the 
curves corresponding to w/k B T = 2 in the insets of Figs. 8 (a) and 8 (b). 

3.2 Approximate solution in the 2D case 

Because the structure of lattice space plays such a fundamental role in deter- 
mining the statistics of /c-mers, it is of interest and of value to inquire how 
a specific lattice structure influences the main thermodynamic properties of 
adsorbed /c-mers. Following this line of thought, we use in this section the 
lattice-gas language again and assume the same model as in Subsection 3.1 
with one exception: the sites form a 2D square lattice instead of a ID lattice. 

In Fig. 9 simulated isotherms are compared to theoretical ones from Eqs. 
(24), (26) and (31) for dimers adsorbed on 2D lattices with different values of 
c and w/kBT: (a), c = 1 and w/ksT = —1,-0.5,0,0.5,1,2; (b) c = 10 and 
w/k B T = -2, -1, 0, 1, 2; and (c) c = 100 and w/k B T = -2, -1, 0, 1, 2. 

In the attractive case, the two theoretical approximations agree qualitatively 
well and the adsorption isotherms for BWA (dashed lines) and QCA (solid 
lines) are hardly distinguishable from each other. The differences between nu- 
merical and theoretical results can be much easily rationalized with the help of 
the absolute error, A9(p/p Q ), which is defined as A9{p/p ) = \9 t heor — 9 S im\p/ Po /9, 
where 9 sim (9 t h eor ) represents the surface coverage obtained by using MC sim- 
ulation (analytical approach). Each pair of values (9 sim , 9 theor ) is obtained at 
fixed p/po- The curves of A9(p/p Q ) vs p/p (data are not shown here for sake of 
simplicity) indicate that, in all cases, QCA leads to appreciably better results 
than BWA. 

Note that the stronger the lateral interaction, the more steep the adsorption 
isotherm becomes. This behavior could be indicative of the existence of a first- 
order phase transition at low temperatures. However, the study of the critical 
behavior of the system is out of the scope of the present work and will be 
object of future studies. 

With respect to repulsive interactions, the differences between QCA and BWA 
are very appreciable. Beyond quantitative discrepancies, there exists qualita- 
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tive differences between both approximations. Thus, while QCA is practically 
independent of c, the discrepancies between BWA and MC results diminish 
appreciably for larger values of c. 

Summarizing, QCA gives a much better description of the 2D MC adsorption 
isotherms than the BWA. In the particular case of repulsive interactions, the 
disagreement between MC and BWA turns out to be significantly large, while 
QCA appears as the simplest approximation capable to take into account the 
main features of the multisite-occupancy adsorption. In fact, there exists a 
wide range of w/k B T's (—2 < w/k B T < 2), where QCA provides an excellent 
fitting of the simulation data. In addition, most of the experiments in surface 
science are carried out in this range of interaction energy. Then, QCA not 
only represents a qualitative advance in the description of the multilayer ad- 
sorption of /c-mers with respect to the BWA, but also gives a framework and 
compact equations to consistently interpret thermodynamic adsorption exper- 
iments of polyatomics species such as alkanes, alkenes, and other hydrocarbons 
on regular surfaces. 

As indicated in the previous section, it is of interest to study the behavior 
of 2D multilayer isotherms (with k > 2) in the low-pressure region in com- 
parison with BET isotherm. For this purpose, the adsorption isotherms are 
plotted in the low-pressure regime and fitted with the linearized form of the 
BET equation. A typical example is shown in Fig. 10. The values of the pa- 
rameters used in the figure are: k — 2, c = 10 and w/ksT = —1,0, 1. Open 
symbols represent Monte Carlo data, full symbols correspond to theoretical 
results obtained from QCA and lines correspond to linearized forms of the 
BET equation. As discussed in Fig. 9, QCA provides very good results in the 
limit of attractive lateral interactions and its accuracy diminishes for repulsive 
ad-ad interactions. 

The analysis of Fig. 10 was repeated for different values of c and w/ksT. In all 
cases, the monolayer volume was calculated from the slope and intercept of the 
linearized form of the BET equation. The results are shown in Fig. 11. Several 
conclusions can be drawn from the figure: (1) QCA agrees very well with 
the numerical results in all range of c and w/k B T studied; (2) the differences 
between v m and vbet diminish as the parameter c increases; and (3) attractive 
lateral interactions compensate the effect of the multisite occupancy. In other 
words, the behavior of vbet/v„i vs w/ksT is similar to the one described 
above for the ID case. 

Finally, the effect of the k-mer size on vbet / v m is analyzed in Fig. 12. In the 
figure, simulation results for k = 2 (full symbols) are compared with the corre- 
sponding ones obtained for k = 4 (crossed symbols). As can be observed, the 
behavior of vbet/v„i vs w/ksT does not significantly vary as the /c-mer size 
changes from k — 2 to k — 4. In addition, the agreement between QCA and 
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MC data (not shown here for sake of clarity) remains very good. Even though 
MC simulations of larger linear adsorbates on regular 2D lattices would be nec- 
essary to confirm the applicability of Eqs. (24) and (31), it should be pointed 
out that QCA is a good analytical approach considering the complexity of the 
physical situation which is intended to be described. 



4 Conclusions 



In this work, we have studied the multilayer adsorption of interacting poly- 
atomic molecules. Two analytic isotherms were developed in the framework 
of the BWA and the QCA. The polyatomic character of the adsorbate was 
modelled by a lattice-gas of /c-mers. With respect to lateral interactions, the 
ad-ad couplings in the monolayer were explicitly considered in the solutions. 
The range of validity of both isotherms was analyzed by comparing theoretical 
and MC simulation results. 

The new formalism from QCA leads to exact results in ID and provides a 
close approximation to study multilayer adsorption of interacting polyatomics 
on 2D surfaces. On the other hand, the artificial effects that the BWA induces 
on the main thermodynamic functions can now be rationalized and compared 
with other analytical approaches. In this sense, we have shown that the dis- 
agreement between BWA and MC simulations increases as (i) the temperature 
is decreased (or the ratio w/ksT is increased) and (ii) the /c-mer size is in- 
creased. 

In addition, we have studied the ID and 2D BET plots obtained using the 
analytic and simulation isotherms. For non-interacting /c-mers, we found that 
the use of BET equation leads to an underestimate of the true monolayer 
volume: this volume diminishes as k is increased. The situation is different 
for the case of interacting molecules. Thus, attractive lateral interactions fa- 
vor the formation of the monolayer and, consequently, compensate the effect 
of the multisite occupancy. In this case, the monolayer volume predicted by 
BET equation agrees very well with the corresponding true value. In the case of 
repulsive couplings, the lateral interactions impede the formation of the mono- 
layer and the BET predictions are bad (even worse than those obtained in the 
non- interacting case). Both the compensation effect for attractive interactions 
and the underestimation of the monolayer volume for repulsive interactions 
are more important for 2D systems. 
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6 Appendix 1 



In order to determine c n , we will use a similar scheme to that employed in 
Appendix A of Ref.[18]. Here, we restrict the analysis to ID systems. 

We start calculating the inflection point of the adsorption isotherm 

d 2 e 



d 2 (p/po) 



and p/p — > (low density). (37) 



By calculating the second derivative of 9 in Eq. (24), we obtain: 

' "(i-p/po) + (i-p/ P o) 2 (i-p/por 1 ] 



where 9" and 9' represent 89 / d{p/p ) and d 2 9/d(p/p ) 2 , respectively. Now, by 
taking lim p / po ^o m Eq. (38), which implies lim^^o, the following relation is 
obtained: 

= 6l + 26' 1 . (39) 



The last equation allows us to obtain c n from the monolayer adsorption isotherm. 
At low density (#i — > 0), the virial expansion approach can be considered as 
an exact result [29]. As usual, Ai(#i) can be written as a power series. Thus, 
by using Eq. (23), 

Ai + c 2^i=o + c 



Differentiating both sides of the last equation with respect to p/po, we obtain 

= \' 1 (\ 1 + c) + \ 1 \' 1 \[c 

( Al + C )2 (Ai + c) 2 1 J 

(E^o^i + c) 2 



Note that Ai — > as 9i — > and, consequently, a = 0. By taking the limit as 
0i -> 0, Eq. (42) results 

1 = ^ (43) 
(c„) 
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and 

0[ 



CJi 



Now, by calculating the second derivative of Eq. (40), we obtain 

_ cAgAi + c) 2 - cA|2(Ai + c) _ c(Ai + c) , , 2 x 

" (aTT^ " w+w { 1 1 1 " 1 J 

A'/c 2 - 2A; 2 c A'/c - 2A/ 2 



= 



c 3 c 2 



= A'/c - 2\[ 2 
being 

oo 



i=0 



and 



i=0 



By taking the limit as 6i — > 0, Eqs. (48) and (49) can be written as: 
Ai(0x - 0) = ai0j 
A^^ -> 0) = 2a 2 0f + a^i'. 



Then, by introducing Eqs. (50) and (51) in Eq. (47) and by using Eq. 
we obtain 

2 2 

= -2a 2 % + c n a x e'[ + 2c„a 2 % 
af af 

c 2 

= -2c 2 + c n ai0" + 2c n a 2 ^, 



and 



o 2 _ o Q2C n o 2 
^1 ~ ~ - Z3~- 



c n a\ a\ a 
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Finally, by introducing Eqs. (44) and (54) in Eq. (39), c n can be written in 
terms of a\ and a 2 : 







2c n 2a 2 cl , 2c„ 



ai af 



+ 



ai 



(55) 



and 



2a? 
a 2 ' 



(56) 



Now, we start with the calculation of ai and a 2 [29]. For this purpose, we write 
the first two terms of the grand partition function of a lattice-gas of N k-mers 
on M sites, 



E k (M, Ai) = 1 + Q k (M, l)Ai + Q fc (M, 2)A? 



(57) 



where Q k (M, 1) and Qk(M,2) represent the partition functions for one and 
two k-mers, respectively. 



By using the following expression, 



In [1 + /(:r)] = In [1 + /(0)] + 



/'(*) 



1 + /(x) 



x=0 



1 f"(x) f"(x)f(x 



-3- 



3! 1 + f(x) [1 + f( x )] 



+ 2 



1 / /" 


» 




[ /'(*) 


2! { 1 + f(x) 


Li+/(x) 


\ /'(*) - 




x 3 + • • • , 

z=0 


Li +/(*)_ 





where /(x) is an arbitrary function, f'(x) = df/dx, f"(x) = d 2 f/dx 2 and 
f"'{x) = d 3 f/dx 3 , lnS fc (M, Ai) can be written as 



In E k (M, Ai) = Q fc (M, l)Ai + - [2Q fc (M, 2) - Q 2 k (M, 1)J A? H (59) 



and the surface coverage results 



M M 9Ai 



M,fe 



= ^Qk(M, l)Ax + A [2Q fc (M, 2) - g^ (M; i) 



M 

Mi + 6 2 Ai + 



z=0 



(60) 
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where the frj's are the well-known virial coefficients. In this case, 
h = ^Q k (M,l) 



(61) 



and 



k 



h = — [2Q k (M,2)-Q k (M,lf 



(62) 



In addition, the relationship between the V s and the a^s can be obtained by 
simple algebra. Thus, 



ai = 



(63) 



and 



a 2 = -To 



b_2 



(64) 



On the other hand, Q k (M, 1) and Q k (M, 2) can be calculated as 

Q k (M, 1) = M (65) 

and 



g fc (M,2)=^ e- 0/fcsT + ^e-^ T 



(66) 



where the configurational factors can be easily obtained. Thus, 



g — —M [M — {2k + 1)] 



q x = —M 2. 
2! 



(67) 
(68) 



Then, 



k 



h = ^ k (M, 1) = k, 



1 1 



(69) 
(70) 
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b 2 = ^[2Q k (M,2)-Q k (M,l) 

= {M [M - (2k + 1)] + 2Me~ w/kBT - M 2 } 
= 2ke~ w/kBT - (2k + l)k, 

2k + l- 2e- w / kBT 



a 2 = — , 3 



k 2 



(71) 
(72) 



and, finally, 



2a? 
a 2 



2k + 1- 2e- w l k B T 



(73) 



In the case w — 0, Eq. (73) reduces to the expression obtained previously for 
noninteracting admolecules [18]. 



24 



References 

[I] T. L. Hill, An Introduction to Statistical Thermodynamics, Addison Wesley 
Publishing Company, Reading, MA, 1960. 

[2] A. Clark, The Theory of Adsorption and Catalysis, Academic Press, New York 
and London, 1970. 

[3] W. A. Steele, The interaction of gases with solid surfaces, Pergamon Press, New 
York, 1974. 

[4] A.W. Adamson, Physical Chemistry of Surfaces, John Wiley and Sons, New 
York, 1990. 

[5] W. Rudzihski and D. Everett, Adsorption of Gases on Heterogeneous Surfaces, 
Academic Press, New York, 1992. 

[6] I. Langmuir, J. Am. Chem. Soc. 40 (1918) 1361. 

[7] S. Brunauer, PH. Emmet, E. Teller, J. Am. Chem. Soc. 60 (1938) 309. 

[8] S.J. Gregg, K.S.W. Sing, Adsorption, Surface Area, and Porosity, Academic 
Press, New York, 1991. 

[9] P. J. Flory, J. Chem. Phys. 10 (1942) 51. P. J. Flory, Principles of Polymers 
Chemistry, Cornell University Press, Ithaca, NY, 1953. 

[10] M. L. Huggins, J. Phys. Chem. 46 (1942) 151. M. L. Huggins, Ann. NY. Acad. 
Sci. 41 (1942) 1. M. L. Huggins, J. Am. Chem. Soc. 64 (1942) 1712. 

[II] E. A. Guggenheim, Proc. R. Soc. London A183 (1944) 203. 
[12] E. A. DiMarzio, J. Chem. Phys. 35 (1961) 658. 

[13] T. Nitta, M. Kuro-oka, T. Katayama, J. Chem. Eng. Jpn. 17 (1984) 45. 

[14] W. Rudzihski, K. Nieszporek, J. M. Cases, L. I. Michot, F. Villieras, Langmuir 
12 (1996) 170. 

[15] G. L. Aranovich, M. D. Donohue, J. Colloid Interface Sci. 175 (1995) 492. 
[16] G. L. Aranovich, M. D. Donohue, J. Colloid Interface Sci. 189 (1997) 101. 
[17] J. L. Riccardo, A. J. Ramirez-Pastor, F. Roma, Langmuir 18 (2002) 2130. 
[18] F. Roma, A. J. Ramirez-Pastor, J. L. Riccardo, Surf. Sci. 583 (2005) 213. 
[19] E. Ising, Z. Physik 31 (1925) 253. 

[20] C. Domb, in: C. Domb, M.S. Green (Eds.), Phase transitions and critical 
phenomena, Academic Press, London-New York, 1974, Vol. 3, pp. 1-95; M.E. 
Fisher, Rep. Prog. Phys. 30 (1967) 731. 

[21] L. Onsager, Physical Review 65 (1944) 117. 



25 



[22] H. Bethe, Proc. R. Soc. London A150 (1935) 552. 
[23] P.J. Flory, J. Chem. Phys. 10 (1942) 51. 

[24] P.J. Flory, Principles of Polymers Chemistry, Cornell University Press, Ithaca 
N.Y., (1953). 

[25] A. J. Ramirez-Pastor, T. P. Eggarter, V. D. Pereyra, J. L. Riccardo, Phys. Rev. 
B 59 (1999) 11027. 

[26] A. J. Ramirez-Pastor, J. L. Riccardo, V. Pereyra, Langmuir 16 (2000) 10167. 

[27] M. Davila, F. Roma, J. L. Riccardo, A. J. Ramirez-Pastor, Surf. Sci. 600 (2006) 
2011. 

[28] N. Metropolis, A.W. Rosenbluth, M.N. Rosenbluth, A.H. Teller, E. Teller, J. 
Chem. Phys. 21 (1953) 1087. 

[29] F. Roma, A. J. Ramirez-Pastor, J. L. Riccardo, Langmuir 19 (2003) 6770. 



26 



Table 1 

Values oiwrninlksT and n (see discussion of Figs. 4 and 5) for different k-mer sizes. 



k-mer size 


Wrnin/ kjjT 




k = 1 


S3 1.1 


» 2.23 


k = 2 


« 1.75 


« 3.31 


fc = 4 


~ 2.5 


» 5.24 


fc = 10 


» 3.5 


~ 11.26 
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Fig. 1. Adsorption isotherms for fe-mers on ID lattices, c = 1 and different values 
of w/ksT (as indicated), (a) k = 1; (b) k = 2 and (c) k = 4. Symbols, solid lines 
and dashed lines represent results from Monte Carlo simulations, QCA and BWA, 
respectively. 
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Fig. 2. As Fig. 1 for c = 10. A zoom of the low-pressure region is presented in the 
inset of each figure. 
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Fig. 3. As Fig. 1 for c = 100. A zoom of the low-pressure region is presented in 
inset of each figure. 
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Fig4a: Garcia el ai. 




1.25 




p Fig4c: Garcia el al. 



Fig. 4. (a) Coordinates of the point of inflection (being Xp and Yp coverage 
and relative pressure, respectively) for monomers (k = 1) adsorbed on ID 
lattices and different values of w/k B T (w/k B T = 2,1.5,1.0,0.5,0.25,0.0,-0.05, 
-0.25,-0.3,-0.4,-0.45,-0.55,-0.75,-0.95,-1.1,-1.2). As a reference, the 
curve corresponding to wjksT = is highlighted. Each point on a given curve 
corresponds to a determined value of c. Solid circles indicate the values of c 
where the inflection point disappears, (b) Same as in part (a) for dimers (k = 2) 
adsorbed on ID lattices and different values of w/ksT {w/ksT = 2,1.5,1.0,0.5, 
0.0, -0.1, -0.25, -0.4, -0.45, -0.5, -0.5$-0.6, -0.65, -0.7, 

-0.75,-0.8,-0.9,-1.0,-1.25,-1.5,-2.0). (c) Same as in part (a) for 
tetramers (k = 4) adsorbed on ID lattices and different values of w/ksT 




Fig. 5. lnc m (as it is indicated in the text) as a function of w/UbT for different 
values of k (k = 1, 2, 4 and 10) and ID lattices. From the slope of the curves in the 
range w/ksT > one obtains k (see discussion in the text). 
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Fig. 6. ( function of w/ksT for k = 1. The meaning of the solid lines and 

the symbols is explained in the text, (b) As part (a) for k = 2 and (c) As part (a) 
for k = 4. 
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0.0-| , 1 , 1 , 

0.0 0.1 0.2 0.3 

P/P 



Fig. 7. \vlVo\l\p (1 —p/Po)] versus p/p for a typical case (c = 10 and k = 2) 
and different values of w/ksT (as indicated). All curves are plotted in the range 
(0 — 0.3) of relative pressure and v m is set equal 1 (in arbitrary units). Inset: 
[p/Po]/[ v (1 — p/Po)] versus p/p for k = 1, c = 10 and three different values 
of w/ksT: full circles, w/ksT = 0; open squares, w/ksT = —1 and full squares, 
w/k B T = 1. 
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Fig. 8. (a) Dependence on w/ksT of the monolayer volume obtained by using BET 
analysis for multilayer k-mer adsorption on ID lattices and c = 10. The curves 
correspond to different values of k as indicated. In the inset, the data are plotted 
as a function of k for three values of w/ksT (w/ksT = —2, w/ksT = and 
w/ksT = 2). (b) Same as part (a) for c = 100. 
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Fig. 9. Comparison between theoretical and simulated adsorption isotherms for 
dimers adsorbed on square surfaces and different values of w/ksT (as indicated), 
(a) c = 1; (b) c = 10 and (c) c = 100. Symbols, solid lines and dashed lines represent 
results from Monte Carlo simulations, QCA and BWA, respectively. 
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Fig. 10. Comparison between theoretical and simulated adsorption isotherms for 
dimers adsorbed on square surfaces with k = 2, c = 10 and different values of 
w/ksT as indicated. The isotherms are plotted in the range of low-relative pressure. 




Fig. 11. Dependence on w/ksT of the monolayer volume obtained by using BET 
analysis for multilayer A;-mer adsorption on square lattices and k = 2. Open sym- 
bols represent Monte Carlo data and full symbols correspond to theoretical results 
obtained from QCA. The curves correspond to different values of c as indicated. 
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Fig. 12. Dependence on w/ksT of the monolayer volume obtained by using BET 
analysis for multilayer fc-mer adsorption on square lattices and two different values 
of k: k = 2 (full symbols) and k = 4 (crossed symbols). The curves correspond to 
different values of c as indicated. 
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